Feshbach Resonance in Optical Lattices and the Quantum Ising Model 
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Motivated by experiments on heteronuclear Feshbach resonances in Bose mixtures, we investigate 
s-wave pairing of two species of bosons in an optical lattice. The zero temperature phase diagram 
supports a rich array of superfluid and Mott phases and a network of quantum critical points. This 
topology reveals an underlying structure that is succinctly captured by a two-component Landau 
theory. Within the second Mott lobe we establish a quantum phase transition described by the 
paradigmatic longitudinal and transverse field Ising model. This is confirmed by exact diagonaliza- 
tion of the ID bosonic Hamiltonian. We also find this transition in the homonuclear case. 

PACS numbers: 67.85.Hj, 67.60.Bc, 67.85.Fg 



Introduction. — With the advent of ultra-cold atomic 
gases, the fermionic BEC-BCS crossover has been the fo- 
cus of intense scrutiny 0-0]. Tremendous experimental 
control has been achieved through the use of Feshbach 
resonances, which allow one to manipulate atomic inter- 
actions by a magnetic field. By sweeping the strength of 
attraction, one may interpolate between a BEC of tightly 
bound molecules, and a BCS state of loosely associated 
pairs. More recently, the analogous problem for a sin- 
gle species of boson has been the subject of theoretical 
investigation jiJ-Q. An important distinction from the 
fermionic case is that the carriers themselves may con- 
dense. This leads to the possibility of novel phases and 
phase transitions, with no fermionic counterpart. 

In parallel, there has been a significant experimental 
drive to study Feshbach resonances in binary mixtures of 
different atomic species. An important catalyst is the 
quest for heteronuclear molecules as a route to dipo- 
lar interactions and exotic condensates [HI. Recently, 
heteronuclear molecules have been created in 85 Rb- 87 Rb 
bosonic mixtures through both s-wave 



The Model. — Let us consider a two-component Bose 
gas with a "spin" index i, f . We assume that these com- 
ponents may form molecules, m. The Hamiltonian 



ia (ij) a 



h.c. 



E 



(i) 



■Hp 



describes bosons, «,•„, hopping on a lattice with sites i, 
where a =i, f, m labels the carrier. Here, e Q are onsite 
potentials, t a are hopping parameters, (ij) denotes sum- 
mation over nearest neighbor bonds, and U aa > are inter- 
actions. We assume that molecule formation is described 
by the s-wave interspecies Feshbach resonance term 



Hf =9 



h.c. 



(2) 



resonances 



12, m 



Similarly, bosonic 41 K- 87 Rb mix- 
tures have been studied in harmonic potentials (Lj | and 
in optical lattices 1151 . and interspecies resonances have 
also been achieved [16Hl8| . The enhanced longevity of 
molecules in optical lattices [l|3], and the sympathetic 
cooling of one species by the other, make these attractive 
for experiment. Multiple species also provide additional 
"isospin" degrees of freedom, and offer exciting possibil- 
ities for interesting phases and quantum magnetism [20| . 

Motivated by these significant developments we con- 
sider the s-wave heteronuclear Feshbach problem for two- 
component bosons in an optical lattice. Our primary 
goal is to establish and explore the rich phase diagram, 
which supports distinct atomic and molecular supcrfiuids 
in proximity to Mott phases. We also shed light on the 
nature of the quantum phase transitions, and a key find- 
ing is a transition described by the paradigmatic quan- 
tum Ising model. This ubiquitous model plays a central 
role in a variety of quantum many body contexts, and 
a controllable realization in cold atomic gases may open 
new directions on dynamics and frustrated lattices. 



Hi and p-wave species, and : m 



For recent work on the p-wave problem see Ref. 21 1. 
Normal ordering yields : n,i a ni a := ni a (rii a — 1) for like 
>.' '■= Tii a ni a i for distinct species. To 
aid numerical simulations we consider hardcore atoms 
and molecules and set U m ^ — U m ± = U and U-^i = V. 
We work in the grand canonical ensemble with = 
H — [It — Nd , where N T = X^( 71 Jt+ n U + 2n mi) and 
Nn = X)i( n 'T ~ n ii) commute with H . These are the to- 
tal atom number (including a factor of two for molecules) 
and the up-down population imbalance respectively. 

Phase Diagram. — To establish the phase diagram it is 
convenient to first examine the zero hopping limit. This 
helps constrain the overall topology and provides some 
orientation for the general problem. With three species 
of hardcore bosons we need to consider eight possible 
states in the occupation basis, |n^,n T ;n m ). The Fesh- 
bach coupling, g, only mixes |1,1;0) and |0,0;1), and 
the resulting eigenstates |±) have energies E± = e m — 
h/2± a/(/i/2) 2 + g 2 , where h = e m -e±-ef-V and the 
chemical potentials are absorbed into = e-f — /it — Pd , 
ej. = ei - /Lt T + /Lt D , e m = e m - 2/x T . The remaining six 
states have energies E(n^, n^; n m ) = J2 a ^q ?1 q + V n t n i + 
Un m (n^ + nj,). Minimizing over all states we obtain the 
zero hopping diagram shown in Fig. [1] In analogy with 
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FIG. 1: Zero hopping phase diagram showing eigenstates in 
the occupation basis \ni,nf,n m ) with ej, = e-f = {/ = 1, 
1/ = 1.5 and (a) # = 1, e m = 3.5, (b) g = 1, e m = 5, (c) 
g = 2, e m = 3.5, (d) 5 = 1, e m = 2. 



the single band Bose-Hubbard model [22J the total den- 
sity pt = n-i + n-f + 2n rn is pinned to integer values and 
increases with /it- Increasing (decreasing) //d favors up 
(down) atoms as is evident from the definition of Afo • In 
contrast to the single band case [22[ , the topology of the 
diagram changes as a function of the parameters. In par- 
ticular, the shape and extent of the central Mott lobe |— ) 
depends on the strength of the Feshbach coupling, g, and 
may terminate directly with the vacuum state and/or the 
completely filled state as shown in Fig. D](b)-(d). 

We now consider the effect of the hopping terms. To 
decouple these we make the mean field ansatz (f> a = (cii a ) 
for each component, and replace a,i a — > <j> a + (<Zj Q — <f) a ). 
This yields the effective single site Hamiltonian 



(a'lfia + (j)* a a a - \(f> a \' 



(3) 



where Hq is the single site zero hopping contribution to 
([1]), and z is the coordination. We minimize Q to obtain 
the phase diagram in Fig. [5J where the symmetry under 
/it — > 2 — /*t reflects invariance of the Hamiltonian upon 
particle-hole and interchange operations, a a f-> aj,, /zt — > 
f-7-n + U — /*t, (14, <-> Of, when t± = and h = 0. 

The phase diagram has a rich structure and exhibits 
single component atomic and molecular condensates, and 
a region with all three supcrfluid. A notable absence is a 
phase where just two components are supcrfluid. This is 
a consequence of the Feshbach term (0); condensation of 
any two variables acts like an effective field on the remain- 
ing species and induces three component superfluidity. 
In contrast, condensation of a single variable no longer 
acts like a field and single component superfluids are sup- 
ported. For example, with /id > we have ej. > e-f- and 
on leaving the vacuum we enter either the up or molec- 
ular supcrfluid. The former is favored at large hopping 
due to the chosen hopping asymmetry; see Fig. [21 

The appearance of single component atomic superflu- 
ids, with (ai) 7^ or Om) 7^ 0, distinguishes this from 
the homonuclear case [4|, [5[. In the latter, single com- 
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FIG. 2: Center: Mean field phase diagram for /in = 0.4, 
e t = ex. = 1, e m = 2, g = 1, U = V = 0, t t = U = t, 
and t m = t/2. We indicate the one-component up, down and 
molecular superfluids, by SFf, SF^ and SF m , while MI de- 
notes the Mott insulating phase. Phase 3SF has all three com- 
ponents superfluid. We denote first order (continuous) tran- 
sitions by double (single) lines. Junctions between phases are 
indicated by a dot, and the termination of first order lines by a 
cross. Top: Magnified portion of the central lobe tip showing 
the retreat of the first order transitions and the emergence of 
a tetracritical point. Bottom: Magnified portion of the lower 
left region showing the underlying tetracritical points as fin 
is varied. 



ponent atomic superfluids arc absent due to the reduced 
form of the Feshbach term ~ g{a\ n aa + h.c). This dif- 
ference also shows up in the symmetry classification of 
the phase transitions. In the heteronuclear problem de- 
fined by equation ([2]), molecular condensation leaves a 
U(l) symmetry intact (a^ — > e^a^, — > e~ la a,i£) in 
contrast to the Z2 symmetry (a — > —a) of Refs. [J, IE]. 
The transition from the molecular to three component 
supcrfluid is thus expected to be in the XY universality 
class along its continuous segment rather than Ising. 

Landau Theory. — The phase diagram in Fig. [2] dis- 
plays an elaborate network of quantum critical points and 
phase transitions. This topology reveals an underlying 
structure that is succinctly captured by Landau theory. 
In the absence of competition from other phases the lo- 
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cus of continuous transitions from the Mott states to the 
one-component supcrfluids may be determined analyti- 
cally. This may be done using second order perturbation 
theory to locate the \(t> a \ 2 terms, or through diagonaliza- 
tion of when two (f> a are set to zero; e.g. the transition 
from |0, 0; 0) to the up-superfluid with (a^) = <fa ^ oc- 
curs along a segment of ztf = e-j-. More generally, 



E = E ° + \H + \ (CA^t + h.C.) 

a 

+ \j2 x "p\ ( t > a \ 2 \M 2 + o(<t> (i ), 



(4) 



where X aa > 0, and the detailed form of the coefficients 
(but not the structure) depend on the unperturbed Mott 
state with energy Eq. This is a Bosc-Hubbard U(l) 
Landau theory for each component (22j . supplemented 
by couplings allowed by the U(l) x U(l) symmetry of 
(TTJ). A similar model (without permissible biquadratic 
terms) was put forward on phcnomcnological grounds to 
describe pairing in the two-band Bose-Hubbard model 
without a Feshbach term [23| . One may gain insight into 
the Landau theory (j4]) by reduction. Near the lower 
tctracritical point in Fig. [2] for example, <fi^ is massive 
(mj, > 0) and may be eliminated using its saddle point: 

£ = £o + ^]Tm a |</g 2 + ^A Q ^ Q | 2 |0 (3 | 2 + ..., (5) 
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FIG. 3: (a) A pair of atoms, aj,, at, is represented by spin 
down |4].) and a molecule, a m , by spin up |-f|-). Hopping favors 
Neel order, whilst the Feshbach term acts like a transverse 
field, (b) mN 1 ^ 8 versus g obtained from exact diagonalization 
of the ID bosonic model |l} for different system sizes, N at 
density pt = 2. We take = e-f = 1 and e m = 2 (h = 0) and 
U = 1, V = 0. We set t = t t = t v = 2t m = 0.01. The crossing 
point, g c ~ 1.041 x 10 -4 , lies slightly below the critical value, 
g c = J/4, of the Ising model ([6]). (c) Scaling collapse as a 
function of (g — g c )N with the critical exponents ft = 1/8 and 
v = 1 of the 2D classical Ising model. 



where a, /3 =f, m, A aa = X aa and A t ,„ = A tm - 7 2 / m i- 
The behavior of this reduced Landau theory ([5]) depends 
on the sign and magnitude of A^ m . For A^A mm > A| m 
it yields a tetracritical point 24J , whilst away from this, 
and for A^ m < 0, it yields two tricritical points separated 
by a hrst order transition. This evolution is borne out 
in Fig. [21 where we track the development of the phase 
diagram as a function of /^d- Indeed, the whole phase 
diagram may be understood within this reduced frame- 
work as the evolution of three such tetracritical points 
and their particle-hole reflections by eliminating <fo, 
and 4> rn in turn. 

Mott Phases. — Having surveyed the overall phase di- 
agram we turn our attention to the Mott states. This 
will reveal Ising transitions in both the heteronuclear 
and homonuclear lattice problems, unreported in Rcfs. 
[1, Q . To explore the central lobe where the Feshbach 
term is operative, we adopt a magnetic description. With 
a pair of atoms or a molecule at each site, we intro- 
duce effective spins | JJ.) = |1,1;0), | ff) = |0, 0;1); sec 
Fig. [3] (a). The operators S + = aj^a^-aj., S~ = ajaja TO , 
and S z = (n m — nj.nf) /2, form a representation of su(2) 
on this reduced Hilbert space, and deep within the Mott 
phase we perform a strong coupling t/U expansion [25| : 



where T = 2g and J 



H 



(ij) 



(6) 



v _ v -r 2u) > is antifcr- 
romagnetic exchange. We have omitted the constant, 
e rn — h/2 — Jz/8 per site. This takes the form of a quan- 
tum Ising model in a longitudinal and transverse field. 
The longitudinal field, h, reflects the energetic asymme- 
try between a molecule | f|-) , and a pair of atoms 1 4J-) ■ The 
transverse field, T = 2g, encodes Feshbach conversion; 
see Fig. [3] (a). Since XY exchange involves interchanging 
two atoms and a molecule it enters at t 3 /U 2 and may be 
neglected. 

Numerical Simulations. — The model ([5]) is of consid- 
erable importance in a variety of contexts, and under- 
pins much of our understanding of quantum magnetism 
and quantum phase transitions. To verify this realiza- 
tion in our bosonic model, we perform exact diagonal- 
ization on the ID quantum system ([T]) with periodic 
boundaries. At present numerical techniques are con- 
siderably less advanced for multicomponent bosonic sys- 
tems, and the large Hilbert space cx 2 3N restricts our 
simulations to N = 8 sites. We begin with h — before 
exploring finite fields. Owing to the absence of sponta- 
neous symmetry breaking in finite size systems, the stag- 
gered magnetization vanishes. Instead, it is convenient 
to focus on m = (| E^-l) 1 ^ |)/A [H, where Sf = 
[n-im — {nit + n it)/2] /2. Adopting the finite size scaling 
form, m = N'' 3 / U m [(g - g^N 1 ^] [H, we plot mA 1 / 8 
versus g for different system sizes, N; see Fig. [3] (b). 
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FIG. 4: (a) Total SF fraction, /„ [27j], for N = 8 sites and 
the parameters used in Fig. [3] The contour indicates the 
approximate location of the MI-SF boundary, (b) g versus J 
showing the Ising transition g c (J) in the second lobe of the ID 
bosonic model {1} using exact diagonalization. In the small J 
limit, g c = J/4 (solid line) in accordance with the transverse 
field Ising model, (c) /i-dependence of F c = 2g c at fixed J, 
for the ID bosonic (squares) and quantum Ising models (solid 
line) . Here e m = 2 + h and all other parameters as in Fig. [3] 



we establish a quantum phase transition described by 
the paradigmatic quantum Ising model. Potential experi- 
ments include magnetization distributions (2(| , quantum 
quenches [30L l3lj . and response near quantum critical 
points. Realizing such models on frustrated lattices may 
probe connections to dimer models. Finite lifetime effects 
may also lead to complex coefficients, and reveal analytic 
properties such as the Yang-Lee edge j32|. In the light 
of these findings it would be interesting to revisit the nu- 
merical simulations in Rcfs [8, 9]. In particular, we have 
verified that the low excited states of the bosonic prob- 
lems are described by Ising Hamiltonians deep within the 
Mott phase. Significantly, the lack of gapless excitations 
for g > g c , and the suppression of XY exchange, sug- 
gests the absence of novel super-Mott H 0] behavior in 
this region of the phase diagram. It would be valuable 
to explore this. 
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These cross close to the critical coupling, g c = J/4, of 
the purely transverse field Ising model, and the scaling 
collapse is consistent with the critical exponents /3 = 1/8 
and v = 1 for the 2D classical model. Repeating this 
we may track the transition within the Mott lobe. The 
results in Fig. @] (b) show clear Ising behaviour at small 
hoppings. In addition, the finite size corrections to the 
ground state energy yield the central charge, c ~ 0.51, 
consistent with the exact result, c = 1/2. At larger J 
we continue to identify a transition in the Ising univer- 
sality class, but the locus is modified. Our system sizes 
are insufficient to track this all the way to the supcrfluid 
boundary; see Fig. [4] These issues will be addressed in 
detail elsewhere by DMRG. Having established an Ising 
transition at h = 0, we now consider /i / 0. As shown 
in Fig. 0] (c), the location of the critical point evolves in 
accordance with exact diagonalization of ([5]) and DMRG 
results on the Ising model (28|. The leading quadratic 
correction to the ground state energy, SE^ m ~Q.07h 2 N 
a t 9 = 9c [Hj is also compatible [33|. This confirms an 
Ising transition in our bosonic model for a range of pa- 
rameters without fine tuning. This also occurs in the 
homonuclear problem with pt = n a + 2n m = 2 [H, Q. 
The spins are modified due to the occupation factors (e.g. 
S + = m jf aa/\/2) but a transition remains. 

Conclusions. — We have considered the Fcshbach prob- 
lem for two species of bosons in an optical lattice and 
have obtained both a rich phase diagram and the over- 
arching Landau theory. Within the second Mott lobe 
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